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ABSTRACT 


Experimental results have recently demonstrated that an adaptive open loop control strategy can be 
highly effective in the suppression of unbalance induced vibration on rotors supported in active magnetic 
bearings. This algorithm, however, relies upon a predetermined gain matrix. Typically, this matrix is 
determined by an optimal control formulation resulting in the choice of the pseudo-inverse of the nominal 
influence coefficient matrix as the gain matrix. This solution may result in problems with stability and 
performance robustness since the estimated influence coefficient matrix is not equal to the actual 
influence coefficient matrix. Recently, analysis tools have been developed to examine the robustness of 
this control algorithm with respect to structured uncertainty. Herein, these tools are extended to produce a 
design procedure for determining the adaptive law's gain matrix. The resulting control algorithm has a 
guaranteed convergence rate and steady state performance in spite of the uncertainty in the rotor system. 
Several examples are presented which demonstrate the effectiveness of this approach and its advantages 
over the standard optimal control formulation. 


INTRODUCTION 


The active control of unbalance vibration in rotating machinery using magnetic bearings has 
generated a great deal of interest in the last decade [1-10]. Recently, research in this area has focused on 
the application of adaptive open loop (or feedforward) strategies. This type of control has the advantage of 
not placing any constraint on the design of the stabilizing feedback control loop for the rotor. Thus, the 
transient performance of the machine can be optimized without considering its unbalance response. 

An important issue in the application of the adaptive open loop control (AOLC) to industrial 
machines is the stability and performance robustness of the unbalance control algorithm employed. 
Recently, the authors have demonstrated that the stability and performance robustness of an AOLC 
algorithm can be analyzed using structured singular value methods [11]. In this paper, these methods are 
extended to provide a synthesis procedure for the design of the AOLC gain matrix. Several examples are 
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presented which demonstrate the effectiveness of the design procedure and its advantage over the standard 
approach that uses a gain matrix derived from optimal control theory. 

Mathematical Notation 

The two-norm of a vector v is indicated by the notation ||v . The maximum singular value of a 
matrix P is denoted by o(P) and the spectral norm by p(P). The lower and upper linear fractional 
transformations [12] of P are given the notations T,(P,Q) and f u (P,R) respectively where the matrices Q 
and R are assumed to be appropriately dimensioned. The Redheffer star-product [13] of appropriately 
dimensioned matrices P and Q will be denoted by S(P,<2). The structured singular value [12,14] of a 
matrix P is indicated by the notation p 4 (P) . The symbol S A is used to denote the set of all matrices of a 
defined block structure. 


ADAPTIVE OPEN LOOP CONTROL 


Adaptive open loop control has been shown to provide excellent vibration attenuation over an 
operating speed range, to quickly respond to sudden changes in unbalance, and to be computationally 
simple [10]. 

The concept of adaptive open loop control is quite simple. Synchronous perturbation control 
signals are generated and added to the feedback control signals. The magnitudes and phases of these 
sinusoids are periodically adjusted so as to minimize the rotor unbalance response. In this paper, a 
particular form of adaptive open loop control, referred to as convergent control, will be examined. The 
convergent control algorithm uses a model of the rotor system where vibration is related to the applied 
open loop signals via 


X = TU + X 0 

where X is a n-vector of the complex synchronous Fourier coefficients of the n vibration measurements, U 
is a m-vector of the complex synchronous Fourier coefficients of the m applied open loop signals, X 0 is a 
n-vector of the complex synchronous Fourier coefficients of the uncontrolled vibration, and T is a n x m 
matrix of complex influence coefficients relating the open loop signals to the vibration measurements. 

The influence coefficient matrix is the transfer function matrix of the supported rotor (with feedback 
control) from perturbation forces at the bearings to the displacements at the sensors, evaluated at the rotor 
operating speed Q , 


T = G f ^ d (jQ) m C(y'QI - A) -1 B + D 

where the rotor/bearing system is described by the state space model (A,B,C,D): 




A 

B 

C 

D 
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Since the convergent control algorithm updates the control vector U periodically, the subscript i 
will be used to denote the t'th update. The time between the i'th and t'+l'th update is referred to as cycle i 
of the algorithm. During cycle i, the control U i is applied resulting in the vibration X t . This vibration 
vector is computed during cycle i using the measurements from the position sensors. The vibration vector 
is related to the control vector by 


X^TU^X 0i ( 1 ) 

where T is assumed to be changing slowly and therefore is not subscripted. During cycle i the next update 
of the control vector U M must be computed from the available information (i.e., U i and X : ). Convergent 
control uses the control update or adaptation law 

U M -U t +AX t (2) 

where A is a gain matrix. The standard approach for determining the gain matrix A is through the 
formulation of an optimal control problem. Minimizing the quadratic performance function 

J - 4CAT,.,} ( 3 ) 

where £{} is the expected value operator, results in the optimal gain matrix [10] 

a = -[t‘t]~ 1 t' (4) 

For implementation, the matrix T for a particular operating speed can be estimated either on-line or 
off-line. When on-line estimation is employed, the gain matrix A is based upon a recently measured 
influence coefficient matrix. If the estimator can track the changes in this matrix, the estimate will be 
accurate and the robustness of the algorithm is not an important issue. However, as the authors have 

experimentally demonstrated [10], the use of continuous on-line estimation may result in large synchronous 
response when there is a sudden change in rotor speed or unbalance. 

An off-line estimate may be obtained through either (1) the injection of test forces using the 
bearings, or (2) modeling of the machine's dynamics (rotor, amplifiers, sensors, and feedback controller). 
Either of these estimates, T can be considered to correspond to a nominal state space model of the system 
given as follows 


G f^i s ) 


A 

b‘ 

c 

D 


where the estimate is related to this nominal model via the equation 

f-G^ d {jQ) 


If the influence coefficient matrix estimate is used to compute the gain matrix A, the nominal system 
optimal gain matrix results 
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( 5 ) 


Because of the possible inaccuracy of the off-line determined estimate, the robustness of the algorithm 
must be considered. With sufficient error, the nominal system optimal gain matrix may result in unstable 
adaptation. In this paper, an alternative design method is proposed for the determination of the gain matrix 
A. This method results in a gain matrix which has a worst case performance that is guaranteed to be 
within a specified factor of the optimal performance. 

The best synchronous performance that can be obtained through active control (open loop or 
feedback) as measured by the minimum value of the quadratic performance index denoted J opl is 

j = Y* X 

J opt opt n opt 

where the minimal vibration vector X opt is given by the expression 

*„-[/-r(r r)V]jr„ 

When the estimate used in the adaptation of the control schedule, T, is equal to the actual influence 
coefficient matrix, T, this performance is achieved by the convergent control algorithm (with A given by 
Eqn. 5) in one update. When the estimate is in error, the adaptation process, governed by Eqns. (1) and 
(2), results in the control vector either growing unbounded or converging to a control vector U n that may 
not be equal to the optimal control vector. If the control vector converges to a equilibrium vector U n , then 
the adaptation process is said to be stable. A necessary and sufficient condition for stability is 
p(/ + AT) < 1. A sufficient condition for adaptation process stability is given by the following condition 

[9]: 

o(l + AT)< 1 (6) 

This stability condition requires that the distance of the control vector from the equilibrium control vector, 
|| U i -£/ n ||, decrease with each update. This more conservative condition, exponential convergence, is 

required to obtain an upper bound on worst case performance. Therefore, Eqn. (6) will be used as the 
stability criterion throughout this paper. If the adaptation process is stable, the steady state value of the 
vibration vector is given by 

X n =[l-T{AT)-'A]x opl (7) 


ROBUSTNESS 

In application of the algorithm, only an estimate of the matrix T can be used in determining the 
gain matrix A. If this estimate, T, was obtained through testing, it may differ from the actual T due to 
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inaccuracies in the estimation process, and/or to changes in the machine's dynamics since the estimate was 
obtained. If the estimate was obtained from a computer model of the rotor, it may differ due to modeling 
errors. 


Both changes in the machine's dynamics and modeling errors usually can be represented by a 
structured uncertainty representation. That is, several parameters 6 l ,9 2 ,...9 j ,-6 p of the dynamic model 
(e.g., the effective stiffness or damping of a seal) are different from those that produced the influence 
coefficient estimate T. Furthermore, a structured representation of uncertainty indicates how each of these 
parameter affects the elements of the influence coefficient matrix T. Each parameter 0 i may differ from its 

nominal value 9 j (i.e., the value which produced the estimate f) with this difference bounded as follows 

8 ,- 6 ,**, 

Thus, each parameter remains in some known neighborhood about its nominal value. Each of these 
parameters may be either real or complex, and therefore this neighborhood may be either a real interval or 
a complex "ball" centered on the nominal value. 

If the state space model matrices (A,B,C,D) are affinely dependent on the uncertain parameters 

, then the influence coefficient matrix of the rotor/bearing system can be represented by a 
linear fractional transformation (LFT) of the following form [11] 

T =5„(G(;Q), A,) = G 22 (/Q) + G 21 (yQ)A,[/ -G n {jQ)A s \' G i2 {jQ) 

where the nominal influence coefficient matrix is given by 

Ga(M)- G,^(jQ)-T 

and A f is a block diagonal matrix of the uncertain parameters. The uncertainties considered determine the 
particular block structure. The set of all matrices of a given block structure is denoted by S A . Note that 
through appropriate scaling of the matrices G u (jQ) and G 12 (;Q), the uncertain parameters can all be 
considered to satisfy 


Throughout the remainder of this paper, the notation ( jQ ) will be suppressed. It will be understood that 
G tj represents a transfer function matrix evaluated at the operating speed of the rotor. 


Stability Robustness 


The adaptation process is said to have exponential convergence to steady state control vector U n 
with convergence rate e c if 
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This is achieved for all T given by the family of matrices 


T-r.{G, A.) 

if and only if (S u ) < 1 where 

S„-S(G,V) 


A a G5 a : a(A,)sl 



( 8 ) 


( 9 ) 


and the uncertainty structure is given by 


A = 


A, 


A / 


where A, is a structured block representing the parametric uncertainty, A v GS A , and A f is a full complex 
block, A, GC mxm [11]. 


Performance Robustness 


As discussed previously, an error in the influence coefficient estimates may cause a decrease in 
performance (unbalance response attenuation). Performance is measured using a quadratic performance 
index of the steady state vibration 


\\x.f (io) 

As shown by the authors [11], the steady state performance for an uncertain rotor system given by Eqn. (8) 
is bounded as follows 


where 


J. <k 2 J 


opt 


k = j 3 + vi9ct(a) 




a 


In this equation, a is a free parameter (a > &), and /?, v, and & are given by the expressions 



r 

aA A 


"a 1 

' 


miny: j u a (S(G,W)) < 1, W ■ 

Y >0 

“/ '-i 
J Y 

, A = 

A,_ 

, A f GC nxn > 


( 11 ) 

( 12 ) 


(13a) 


156 




f 

‘ G„ 

0,2 ' 



rA, i 


miny: ^ A 

y >0 

l 

j 

-o 21 

y 

) 

A 

> 

II 

. A /_ 

A i — /^ymxn 

, A / EC 




1 


1- miny 

0<y <1 


: ^ a (S(G,V))<1, 


' o r 


rA s i 


La 

j r 

, A = 

A /_ 

, A, EC mxm > 


(13b) 


(13c) 


As Eqn. (12) indicates, the variable k can be decreased so that it is very close to the value of through 
increasing a. For large a, can be considered a good approximation to k. If this approximation holds, 
then ju a (S(G,IP)) < 1 implies that J n < fi 2 J opr The derivation of Eqn. (12) requires that the stability 

condition, Eqn. (6), be satisfied. 


Convergence Rate and Performance 


Both a convergence rate of e c and an (approximate) performance of J n < fi 2 J opl can be 
"guaranteed" if the following condition holds 


A » a ( S )<1 


S(G,V0 


S{G,W) 


(14) 


where the structure of the uncertainty block A is given by 


A = 


A 1 , 


a 1 „a 2 

mxm 


A\ EC 
A 2 f EC'"'" 


(15) 


where V and W are as defined in Eqns. (9) and (13a). Note that both V and W in Eqn. (14) are affinely 
dependent on the gain matrix A. 

An upper bound on the structured singular value / u is available via the equation 


M s ) 


^ 1 min q : q * cr 

VlPkA 1 

\GlGmGr 


(I +^ 2 PKS©;' - jqQ M )(l +QIY* , q<q 


(16) 


where {'D L ,'D R ,g L ,g M ,g R ] must be chosen from sets of appropriately structured block diagonal matrices 
[12] and q is a positive scalar. The minimization problem given in Eqn. (16) is convex, and therefore the 
upper bound can be determined very quickly and easily [15]. 
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SYNTHESIS OF ROBUST GAIN MATRICES 


The result of Eqns. (14) and (16) gives rise to a design algorithm for adaptive open loop control 
gain matrices. If the upper bound on ^t A (S), Eqn. (16), is minimized through choice of A to be less than 

one, then the convergence and performance robustness specifications will be achieved. The actual 
performance robustness must be checked after this design procedure since the performance specification is 
only approximate. In practice, however, a can be chosen to be very large, so that the design method using 
the approximation produces a gain matrix satisfying the desired performance specification. 


The design problem is as follows: 


minimee q : 

,q 

G lG m Gr 


r 


< 7 S °1 


(/+^) 7 K s ( a K' 


V 


WGmX 1 





) 


q<q 


(17) 


While the individual minimization with respect to \p L ,T> R ,Q L ,Q M ,Q R ,q] and with respect to A are both 
convex, the combined minimization with respect to {A, 2) L ,T> R ,Q L ,Q M R ,q} is not convex. Therefore, 
the minimization of Eqn. (17) is handled via an iteration between the two convex problems. First, 
{'D L ,‘D R ,Q L ,g M ,g R ,q] will be fixed (starting with = I, Q L ,Q M ,Q R = 0) and the minimization is 

carried out with respect to A. Then, A is fixed and the minimization is carried out with respect to 
\p l ,T> R , g L ,g M ,Q R ,q) . This second step is just the computation of the upper bound on £t A (S(A)) given 

in Eqn. (16) with q as the resulting upper bound on jU a (S(A)). This iterative synthesis procedure is 
analogous to that used in D-K iteration for /z-synthesis [15]. The minimization is stopped when the upper 
bound on ju a (S(A)) is less than one. 


EXAMPLE PROBLEMS 


Example # 1 

First, the two mass system, shown in Figure 1, is considered. While this system is very simple, its 
vibration behavior is analogous to that of a rotor supported in magnetic bearings. For this example, the 
spring between the masses has a uncertain stiffness, k v with a nominal value of 1. The values of the other 
parameters are: n\ = 05, = 1, k 2 = 1, and c = 1. The "operating speed" of 1.2 radian/s is first 

considered for this example problem. This is nominally between the first and second critical speeds of the 
system. With a 40% decrease in stiffness, the second natural frequency approaches the operating speed. 
The first natural frequency is relatively insensitive to this variation in stiffness. 

A gain matrix was first designed to meet the following requirements: 

with an uncertainty in the stiffness k t of ± 50%. From these specifications, the values fi - V 1 .25 and 
e c = 0.9 were used in the synthesis procedure. The free parameter a was set at 10 5 so as to make the 
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/3 - K approximation accurate. The procedure was successful in synthesizing a gain matrix which meets 
all of the specifications. 

The time response of the adaptive open loop vibration control was simulated using Eqns. (1) and 
(2). A comparison of the behavior of the system was conducted using the nominal system optimal gain 
matrix, Eqn. (5), and the gain matrix obtained by the synthesis procedure. In both cases, the convergent 
control algorithm was started with U = 0 (no open loop control forces applied). With an actual stiffness 
k { =05, the quadratic performance index time history with the nominal system optimal control and with 
the robust control are shown in Figure 2. Both time histories have been normalized by the minimum 
value of the performance index J . Note that the robust control results in very little degradation in the 

converged value of the performance index in comparison with the optimal value. However, the nominal 
system optimal gain matrix resulted in unstable adaptation. 

To fully examine the performance robustness of the robust gain matrix control, 100 simulations 
were conducted with the actual stiffness varying from 0.5 to 1.5. Each performance index history was 
normalized by the optimal performance. Figure 3 shows the results of these simulations. In each case, the 
adaptive open loop control converges to a steady state satisfying J n s 1.257,,^ . 

Example #2 

A more complex rotordynamic system is considered in the second example problem. Figure 4 
shows a diagram of a rotor considered for a boiler feedpump application. The rotor has a length of 2.54 m, 
a diameter of 57 mm, and a mass of 226 kg. Two magnetic bearings are used with proportional-derivative 
feedback control for stabilization. The rotor also has a fluid film bearing. For the purposes of this 
example, the magnetic bearing feedback control as well as the fluid film bearing will be modeled as a 
stiffness and a damping at each bearing location as shown in the diagram (stiffness k v k 2 ,k $ and damping 
c,,c 2 ,c 5 ). Parameters k 3 and k 4 represent the motor's stiffness which is due to both magnetic and fluid 
dynamic effects. The nominal values of each of these parameters and the uncertainty in each considered in 
this example are shown in Table 1. The two lowest critical speeds of the nominal rotor system are 2600 
and 4970 rev/min. The balancing forces are injected at bearings #1 and #3 (the magnetic bearings) so as to 
reduce the vibration at bearings #1 and #3 and the motor locations. These inputs and outputs are indicated 
in Figure 4. For this example, the operating speed of the rotor is 3000 rev/min. 

A robust gain matrix was synthesized to meet the following specifications: 

J.* 

Figure 5 shows the time histories of the normalized performance index for convergent control with the 
robust and nominal system optimal gain matrices. In these simulations, the errors in the nominal 
stiffnesses were: 


<5, = 0525, & 2 = 0525, d 3 = 0.613, and <5 4 = 0.088 N/pm. 

which is within the uncertainty bounds given for the parameters. The nominal system optimal gain matrix 
resulted in unstable adaptation for this variation. However, the robust gain matrix produced very good 
performance. The robustness of the synthesized gain matrix was also examined through 1000 time 
simulations with the system's uncertain stiffnesses chosen as uniformly distributed random numbers 
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within the ranges given in Table 1 . The results of these simulations are shown in Figure 6. In all cases, the 
robust convergent control was stable and produced a steady state performance that was bounded as 
originally prescribed in the synthesis specification. 


CONCLUSIONS 


The robustness of adaptive open loop control algorithms for suppression of synchronous vibration 
can be significantly improved through the application of a simple synthesis procedure. This procedure can 
be used to design gain matrices which have both stability and steady state performance robustness to real 
and complex structured uncertainties. 

As the example problems demonstrate, the nominal system optimal gain matrix usually employed 
in adaptive open loop control can have a significant performance degradation and even instability when 
variations occur in the system parameters. In contrast, the performance of the robustly-synthesized gain 
matrix degraded significantly less from the optimal. Its worst case performance is known and was 
significantly better than that produced by the nominal system optimal gain matrix. 
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Table 1. Nominal Values and Uncertainties for Example #2 



*1 

55 

1.75 

* 2 

55 

1.75 

* 3 

0 

0.875 

K 

0 

0.875 

k 5 

55 

1.75 

Damping (Ns/mm) 


91.3 

0 

C 2 

91.3 

0 

C 3 

91.3 

0 
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0 5 10 15 20 25 30 

Updates 

Figure 2: History of performance indices for simulation with robust and nominal system optimal gain 
matrices, k , = 05 . 
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Updates 

Figure 3: History of performance index for 100 simulations with robust gain matrix, stiffness varies 
from 0.5 to 1.5. 



Bearing #1 Bearing #2 Motor Bearing #3 

(magnetic) (fluid film) (magnetic) 


Figure 4: Rotor model considered in example #2. 
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Figure 5. History of performance indices for simulation with robust and nominal system optimal gain 
matrices, example #2. 



Updates 


Figure 6: History of performance index for 1000 simulations with robust gain matrix, example #2, 
maximum J < 1.247 , . 

n opt 
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